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ABSTRACT 

We present a detailed study of the star formation history (SFH) of the Tucana dwarf spheroidal 
galaxy. High quality, deep HST/ACS data, collected in the framework of the LCID project, allowed 
us to obtain the deepest color-magnitude diagram to date, reaching the old main sequence turnoff 
(.F814 ~ 29) with good photometric accuracy. Our analysis, based on three different SFH codes, 
shows that Tucana is an old and metal-poor stellar system, which experienced a strong initial burst of 
star formation at a very early epoch (~ 13 Gyr ago) which lasted a maximum of 1 Gyr (sigma value). 
We are not able to unambiguously answer the question of whether most star formation in Tucana 
occurred before or after the end of the reionization era, and we analyze alternative scenarios that 
may explain the transformation of Tucana from a gas-rich galaxy into a dSph. Current measurements 
of its radial velocity do not preclude that Tucana may have crossed the inner regions of the Local 
Group once, and so gas stripping by ram pressure and tides due to a close interaction cannot be 
ruled out. A single pericenter passage would generate insufficient tidal heating to turn an originally 
disky dwarf into a true dSph; however, this possibility would be consistent with the observed residual 
rotation in Tucana. On the other hand, the high star formation rate measured at early times may 
have injected enough energy into the interstellar medium to blow out a significant fraction of the 
initial gas content. Gas that is heated but not blown out would also be more easily stripped via ram 
pressure. We compare the SFH inferred for Tucana with that of Cetus, the other isolated LG dSph 
galaxy in the LCID sample. We show that the formation time of the bulk of star formation in Cetus 
is clearly delayed with respect to that of Tucana. This reinforces the conclusion of Monelli et al. 
(2010) that Cetus formed the vast majority of its stars after the end of the reionization era implying, 
therefore, that small dwarf galaxies are not necessarily strongly affected by reionization, in agreement 
with many state-of-the-art cosmological models. 

Subject headings: Local Group galaxies: individual (Tucana dSph) galaxies: evolution galaxies: pho- 
tometry Galaxy: stellar content 
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Nearby galaxies in the Local Group (LG) present a 
variety of different properties in terms of stellar popu- 
lations, gas and metal content, morphological type, and 
mass. Understanding their differences and similarities 
and how these evolved with time can bring important 
insights into our knowledge of the formation and evolu- 
tion of the LG, and of galaxies in general. The recovery 
of their full star formation history (SFH) plays a key 
role, being among the most powerful techniques to dig 
into the mechanisms that drove the evolution of stellar 
systems. Precise estimates of the epochs of star forma- 
tion and their duration give direct information to com- 
pare the observed properties with predictions. This is 
particularly important for investigating the first stages 
of the life of galaxies, because the properties of the stel- 
lar populations during the first few Gyr can be used to 
trace the impact of a variety of different physical mecha- 
nisms. For example, a delay in the onset of the star for- 
mation, or a truncation of the SFH at a very old epoch 
can possibly be used to infer the role of cosmic reion- 
ization (e.g.. IIkeuchil [l986; Rccs 119861: lEfstathioul 119921: 
Babul fe Reesl Il992t Ichiba fc Nathl Il994t IQuinn et al l 
19961: iThoul fe Weinberg! 119961: iKepner et al.l 119971: 
Barkana k Loeblll999t i Bullock et all 120001: iTassis et ail 
2003t IRicotti fc Gnedinl 120051; IGnedin fc Kravtsovl 120061: 



Okamoto fc Frenkl2009l : ISawala et al.l2010l:lMaverl2010n . 

On the other hand, an episodic rather than continu- 
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ous SFH might trace interactions of dwarf galaxies i n 
their orbit around bigger systems (iMaver et al.l 120011) , 
or a possible cycle related to changes in t he interstellar 
medium linked to supernova explosions dCarraro et al.1 
120011 IStinson et all [20071 : IValcke et al.ll2008D . WhaTis 
crucial in these kinds of analysis is the capability to pre- 
cisely and accurately recover how these systems formed 
stars as a function of time, and to characterize their 
chemical evolution. 

Within this general context, we selected a sample of 
isolated dwarf galaxy members of the LG, with the aim 
of deriving quantitative SFHs over their entire lifetime. 
The main objective of the LCID project (Local Cosmol- 
ogy from Isolated DwarfF^) is to study the details of the 
evolution of nearby isolated galaxies, compare their prop- 
erties and help constrain cosmological and galaxy forma- 
tion and evolution models. A general description of the 
project, and an overview of the results, are presented in 
Gallart et al. (in prep). This paper is focused on the 
study of one the most isolated dSph galaxies in the LG, 
Tucana, with particular emphasis on the comparison its 
proper ties with those of the Cetus dSph (j Monelli et al.l 
I2010bfl . 

Understanding the currently isolated nature of Tu- 
cana and Cetus is challenging in the context of galaxy- 
formation theories. Models have been proposed sug- 
gesting that dSph galaxies originate from small, initially 
gas-rich galaxies that lost their gas through tidal and 
ram-pres sure interact ions with their large host galax- 
ies (see iMaverl [20101 for a review). Before the LCID 
project, however, the existing observations of Cetus and 
Tucana did not allow us to draw definite conclusions on 
the exact nature of these galaxies, particularly regard- 
ing the q uestion of how extend ed in time their SFHs 
were. In Monelli e t al.l ([2010b[ ) and in this paper, we 
show that Cetus and Tucana are similar to the oldest 
Milky Way dSph satellites such as Draco, Ursa Minor, 
Sculptor and Sextans, and therefore they do not follow 
the morphology-density re lation observed amon g Milky 
Way dSph satellites (e.g. Ivan den Berghlll999j ). Their 
mere existence, therefore, imposes new constraints on 
dSph galaxy-formation models. 

The first mention of the Tucana dSph a ppears in the 
Southern Galaxy Catalogue published b y iCorwin et al.1 
(p85h . but it was later "rediscovered" bvlLavervl 0. 990). 
The first CMD was presented in lLaverv fc M ighcll 
(1990), based on Anglo- Australian 3.9m V, I images. 
Tucana appeared as an elongated spheroid (e=0.5), and 
was classified as a dE5. Although their photometry 
was quite shallow, lLaverv fe Mighelll ([19901 ) could give 
the first estimates of the distance ((m — M) < 24.75), 
metal content ([Fe/H] = —1.9) and luminosity (My = 
—9.5). Tucana was definitively recognized as a mem- 
ber of the LG, with the peculiarity of being the first 
isolat ed dSph, not linked with either the Milky Way or 
M31. lLaverv fc Mighelll (|1990| ) also inferred that Tucana 
is a predominantly old system, based on the lack of young 
bright blue objects. 

All these preliminary find ings were later substan tially 
confirmed by deeper data. [Castclla ni et al.l (|1996t ) pre- 
sented V, I photometry reaching V ~ 26 mag, which al- 
lowed them to estimate the distance ((m— M )j = 24.72 ± 

14 http: / /www. iac.es/project/LCID 



0.20) based on the tip of the red giant branch (RGB), and 
a mean metallicity of [Fe/H] = —1.56, obtained compar- 
ing the Tucana RGB wi th the ridge line of Galact i c glob - 
ular clusters (GCs) from lDa Costa fc Armandrofll (ll990). 
Interestingly enough, they found that the color width of 
the RGB was larger than expected on the basis of the 
photometric errors. They interpreted this as a signature 
of a " consistent spread" in the metallicity of the Tucana 
stars. 

A d ifferent conclusion was reached by iSaviane et al.l 
(1996), who presented V, I data collected with the 
EFOSC camera mounted on the 2.2m ESO tele- 
scope. A comparison with t he same database from 
iDa Costa fc ArmandrofH (|1990T ) yielded a mean metal- 
licity [Fe/H]~ —1.80, and "no convincing indication of 
an abundance sprea d." An interesting point raised by 
ISaviane et al.l ([19961 ) concerns the absolute dimension of 
Tucana. Using an exponential profile, they estimated a 
core radius of 166 pc, similar to that of the small Milky 
Way and M31 satellites (Mateo 1998). 

Another very interesting characteristic of Tucana is the 
strong ho rizontal branch morp hology gradient first de- 
tected bv lHarbeck et all (|2001| ). together with a double 
peak in the color distribution of the RGB stars. Based 
on these two hints, and on the likely absence of in- 
termediate age stars (supported by th e non-detection 
of C stars, iBattinelli fc Demersl 120001 ) . lHarbeck et all 
(|2001f l suggested that Tucana experienced a very fast 
sel f-enrichment. This p ossibility was strongly supported 
by lBernard et al.l ((2008) on the basis of the properties of 
a sample of ~400 RR Lyrae stars. They detected two dif- 
ferent sub-populations of RR Lyraes having distinct spa- 
tial distributions and different mean luminosities, and 
following two period-amplitude relations with different 
slopes. This supports a scenario with two different sub- 
populations older than 10 Gyr, the more centrally con- 
centrated one being slightly younger and slightly more 
metal rich than the other. 

The paper is organized as follows. fj2] summarizes the 
observations and data reduction strategy. In sect. f}3] 
we present the CMD, while the derivation of the SFH 
is extensively discussed in Ej4] In fj5] we compare the 
properties of the two isolated dSphs Cetus and Tucana. 
Discussion and conclusions are presented in iJB] and 

2. OBSERVATIONS AND DATA REDUCTION 

The data presented in this pape r were collected w ith 
the ACS camera aboard the HST (|Ford et al.lll998D . as 
part of the project The onset of star formation in the 
universe: constraints from nearby isolated dwarf galaxies 
(PID 10505, PI C. Gallart). 

The 32 orbits allocated to Tucana were executed be- 
tween April 25 and 30, 2006. The observations were 
split into visits, each of them including two orbits. One 
FA75W and one F81AW image were collected during 
each orbit, with exposure times slightly different among 
the two orbits of the same visit: 1,070 s and 957 s in the 
first orbit, 1,090 s and 979 s in the second one. Summa- 
rizing, the total observing time on Tucana was 34,560 s 
in F475W and 30,976 s in F8UW. Note that, despite 
the fact that Tucana is the most distant galaxy in the 
LCID sample, the exposure times adopted are the short- 
est. This choice was forced due to the constraints im- 
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Fig. 1. — Stacked, drizzled color image of the Tucana field. North is up and East is left. The field of view is ~3'4x3'4. The image shows 
a clear gradient in the number of stars when moving from the center to the outskirts. A sizable number of background galaxies is visible 
as well. 



posed by the South Atlantic Anomaly. 

The pointing (a = 22 /l 41 m 48M3, S = -64° 25' 15^7) 
was centered on the galaxy. Figure [T] shows a drizzled, 
stacked i mage. The parallel W FPC2 field already pre- 
sented in iBernard et al.1 (j2009T ) was intended to sample 
the halo of Tucana along its major axis, at the position 
(22 h 40 m 52^92; -64° 24' 28^9). The CMD of this field, 
located 6'from the Tucana center, does not show any con- 
vincing evidence that stars of Tucana are still present at 
this distance, since no RR Lyrae stars were detected. 

2.1. Data reduction 

The details of the data reduction strategy, com- 
mon to all th e LCID galaxies, have been presented in 
iMonelli et al.l (|2010bfl . so here we just recall the main 
points. We performed two parallel and independent pho- 
tometric r educt ions, using the DAOPHOT / ALLFRAME 
dStetsonl Il99l and DOLPHOT packages (|Dolphinl 
2000b). We applied both codes to the original _FLT im- 
ages, working on the two chips separately. In the case of 
DAOPHOT, we modeled individual PSFs for each frame, 
using bright and isolated stars in the field. Typically, 400 
stars and a Moffat function with index f3 — 1.5 were used. 

The DOLPHOT code was applied using the ACS mod- 



ule and following the recommended photometry recipe 
provided in the manual for version 1.0.3. 

Individual photometry catalogs for each frame were 
calibrated. We a pplied the calibration provided by 
ISirianni et al.l ()2005D identically to both sets of photom- 
etry. The final CMD is presented in Fig. [5] and discussed 
in section [3] 



2.2. Completeness tests 

The completeness tests were performed by adding syn- 
thetic stars to the individual images, and repeating the 
photometric analysis identically as for the real data. 

In the case of the DAOPHOT photometry the list of 
injected sta rs was created follow ing the prescriptions in- 
troduced bv lGallart et al.l (119961) . using a syntheti c CMD 
created with IAC-star (jAparicio fc Gallartl I2004D . This 
was built adopting a constant star formation rate be- 
tween and 15 Gyr, and a metallicity distribution fol- 
lowing a flat distribution in the range 0.0001 < Z < 
0.005. Approximately 1,500,000 stars were simulated in 
each of the two chips, in iterations of ~ 50,000 stars. 
In the case of DOLPHOT, we followed the method de- 



scribed in lHoltzman et al. (2006), simulating one star at 
a time and covering the range — 1 < Mp^^W ~ Mg,uw < 
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Fig. 2. — CMD of Tucana and Cetus obtained with DAOPHOT. Note that the color of the bluest HB stars is similar in both galaxies, but 
the number of blue HB stars is significantly higher in Tucana than in Cetus. The RGB bumps are also obvious features in b oth galaxies as 
indica ted by the arrows on the right side of the RGBs. The arrow on the left side of the Tucana RGB mark the AGB clump. (Mone lTi et alj 
2d I (la'. 
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5 mag, — 7 < Mp^jsw < 8 mag. The density of injected 
stars is higher in the color-magnitude range where most 
of the observed stars are found. A total of 140,000 stars 
were simulated. We found that, at the level of the TO, 
the completeness is of the order of 80%. 

2.3. Comparison of the two photometry setss 

We detect a small zero point offset for the brightest 
stars in the m F475W (m F475WyD Ao - m F 4 75W ,DOL = 
+0.019 mag, for m F 4 7 ^w < 25) and in the rnpguw 
(m F8 i4w,DAo ~ rn F81 4w,DOL = -0.008 mag, for 
ttifsuw < 24) bands, and a trend as a function of 
magnitude. As alrea dy noticed in the case of Cetus 
(jMonelli et al.l l2010bO . ' moving toward fainter magni- 
tudes, DAOPHOT tends to measure stars slightly fainter. 
However, since the trend with magnitude is similar in 
the two bands, the color difference shows a small resid- 
ual zero point (~ 0.032 mag), but no dependence on 
the magnitude. As al ready noticed in the case of Cetu s 
(jMonelli et al.ll2010bl) and LGS 3 jHidalgo et al.l [20Toh . 
and discussed in IHoltzman et al.l ([2006) in the case of 
the WFPC2, small differences in the photometry are ex- 
pected when using different packages. Most importantly, 
we will show in section 14.51 that the impact of the differ- 
ences in the photometry on the derived SFH is negligible. 

3. THE CMD 

Fig. H left panel, presents the final CMD for Tu- 
cana from the DAOPHOT photometry, calibrated to the 
VEGAMAG system and shifted to the a bsolute plane. 
The ri ght panel shows the CMD of Cetus (|Monelli et al.l 
l2010bD . The comparison of the two CMDs shows strik- 
ing similarities in the different evolutionary features. An 
in-depth comparison of the CMDs and the properties of 
these two galaxies is presented in ij5] 

The Tucana CMD presented here is the deepest ob- 
tained for this galaxy to date. It spans more than ^8 
mag, from the tip of the RGB down to ssl.5 mag be- 
low the oldest main sequence (MS) turn-off (TO), or 
Mpsuw ~4. The absence of a blue young MS and 
the morphology of the old MSTO, at M F8 i4\y ~ 2.5 
mag, are typical features of an old stellar system, with 
no strong recent or intermediate-age star formation. In 
agreement with what is commonly found in stellar sys- 
tems, the plume of relatively faint objects appearing at 
1.5 < Mfsuw < 3, < (M Fi75W - M FSliW ) < 0.5 is 
most likely a population of blue stragglers. The possi- 
bility of some residual star formation occurring in the 
last 3-4 Gyr is discussed in M4.51 The RGB appears as a 
dominant feature, between -4 < Mp^uw < 2 mag. The 
width of the RGB suggests a limited spread in metal- 
licity. The two over-densities located along the RGB at 
M FSUW ~ -0.9,-0.7 and M F475W - M FSUW ~ +1.4 
and highlighted by t wo arrows are the RGB bumps (see 
iMonelli et all l2010al for a detailed study of the RGB 
bump feature in several LCID galaxies). The other 
clear over-density at Mpsuw ~ —0.8 and M F 4 7 $w — 
Mpsi4w ~ 1-2, and joining the RGB from the blue side, 
can be identified with the AGB bump and the subsequent 
AGB. 

The Tucana HB, as previously noted by Harbeck et al. 
(2001), shows a complex morphology, which can be ap- 
preciated in detail in the present, much deeper data. 




Fig. 3. — Tucana CMD with superimposed ZAHBs of different 
metallicities Z=0.0001 (dotted), 0.0003 (dashed), 0.0006 (solid) 
and Y= 0.245 (red), 0.30 (green), 0.35 (blue) are shown. Only 
the theoretical predictions with normal helium give a good repre- 
sentation of the HB observed in Tucana. 



It is well populated from the red to the blue side, to 
Mf4 7 5w — Mf&iaw ~ 0. A few bluer and fainter objects 
seem to continue the HB sequence, but they merge with 
the BSs candidates and it is not possible to draw firm 
conclusions on their nature on the basis of the present 
data. A possible scenario to explain the strong blue com- 
ponent of the Tucana HB is that this feature is populated 
by a helium-enriched population. Helium-rich stars have 
been invoked to explain the properties of the most mas- 
sive Galactic globular clusters which have been proved 
to host more than one population dPiotto et al.l 120051 : 
iCaloi fc D Antonal 120051: iMilone et al.ll2008l |2010T) . The 
effect of the helium enrichment is to produce, at fixed age 
and for M F 4 7 5w — M F si4w > —0.2, a bluer and brighter 
HB. Fig. [3] shows the Tucana CMD with selected ZA- 
HBs from the BaSTI database superimposed, for metal- 
licity Z=0.0001, 0.0003, 0.0006 and Y=0.245, 0.30, 0.35. 
The figure shows that the HB is nicely bracketed by the 
canonical helium ZAHBs, and a metallicity in the range 
0.0003-0.001. Note that this is in excellent agreement 
with the results of the SFH, suggesting a mean metal- 
licity of Z = 0.0005 (see g3J). The helium-rich ZAHBs 
appear too bright to explain any significant number of 
Tucana HB stars. 

4. THE STAR FORMATION HISTORY OF TUCANA 

The SFH of Tucana was derived, similarly to 
the other galaxies of the LCID sample, using the 
two phot ometry sets, two stell ar evolution libraries 
( BaSTFl. iPietrinferni et~aT] [200l and Padova/GirardO, 
iGirardi et al.l I2000D and three di fferent SF H codes: 
I AC- pop dA paricio fc Hid"aJgol [20091) . M ATCH (jPolphinl 
I2002T ). and COLE dSkillman et al.ll2003D . We performed 
the same kind of ana lysis already presented for Cetu s 
(|Monelli et al.1 12010bl) and LGS3 (|Hidako et al.l[20Tol) . 
In particular, we used the three SFH codes together 
with the DOLPHOT photometry in combination with 
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Fig. 4. — Summary of the results obtained with different techniques to derive the SFH. Top: Comparison of the cumulative mass fraction 
derived from SFH obtained with the MinnlAC/IAC-pop method, using the BaSTI and Padova library and applied to the DAOPHOT and 
DOLPHOT photometry. Middle: Same as in the upper panel, but the comparison shows the results the from three different SFH codes, 
applied to the DOLPHOT photometry in combination with the Padova/Girardi Bottom: ip{t) derived for the three solutions presented in 
the central panel. 



the Padova stellar evolution library. The IAC-pop 
solutions were derived using both photometry sets and 
both libraries. 

In the following sections we summarize the fundamen- 
tals of the three SFH methods applied. A comparison of 
the results is given in ij4.41 while the Tucana SFH will be 
discussed in ij4.51 



4.1. The IAC method 

The IAC method to derive the SFH of re- 
solved stellar systems, based on the IAC-star/IAC- 
pop/Min nlAC algorithms, h a s bee n ex tensively pre- 
sented in lAparicio fc Hidalgol (|2009f) and iHidalgo et al J 
(2010). Due to the similarities between the Tucana and 
Cetus CMD, we applied the same criteria explained for 

15 \http://www. oa-teramo.in af.it/BASTI 

16 http://pleiadi.oapd.inaf.it/ 



Cetus in I Monelli et al.1 (|2010bO . The reader is referred 
to that paper for more details and for the discussion of 
the different choices. In summary, we have adopted the 
following parameters and functions: 

i) Model junctions- The same synthetic CMDs con- 
taining 8 x 10 6 stars have been used as for Cetus. These 
were generated using the IAC-star simulatoiFI. using: 

• the BaSTI and Padova stellar evolution libraries, re- 
spectively; 

• a constant tp(t) between and 15 Gyr, and metallicities 
uniformly distributed in the range 0.0001 < Z < 0.005; 

• a binary fraction^ 0=O A wit h mass ratios q > 0.5; 

• the IMF from iKroupal ()2001f ) (N(m) dm = mT a dm, 
where a — 1.3, 2.3for stars with mass smaller and larger 
than 0.5 M , respectively); 

• the set of bolometric corrections from IBedin et al J 
(f2005h . 

The observational errors were simulated in the synthetic 
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CMD using the obsersin code (|Hidalgo et al.ll2010h . 
it) Simple populations and CMD sampling: 
The initial age and metallicity bins defining the simple 

populations are: 

age: [1.5 2 3 ...11.5 13.5] xlO 9 years 
metallicity: [0.1 0.3 0.5 0.7 1 1.5 2]xl0~ 3 

The parametrization of the CMD was performed us- 
ing the same c hoice o f bundles and box size presented in 
IMonelli et all (l2010bD ( see Figure [7]). Moreover, the same 
multiple sampling of both the observed and the model 
CMDs (" dithering approach") has been app lied, using 
the MinnlAC code (jAparicio fc Hidalgo! [20091 ) . Twenty- 
four solutions with twelve different samplings of the sim- 
ple populations and two samplings of the positions of the 
boxes grid each, were calculated. 

Hi) External parameters: 

The true distance modulus, (m — M)q = 24.77(3, was 
estim ated using the mean magnitude of the RR Lyrae 
stars dBernard et al J 120091) . We assumed the extinction 
from iSchlegel et alT (fl9*98). E(B - V) = 0.03, corre- 
sponding to Aptfsw — 0.122 mag and Apsi4W = 0.061 
mag. To estimate the effect of the uncertainties affect- 
ing the photometry and the transformation to the abso- 
lute magnitude plane (calibration, distance, extinction) 
and possible model zero points, MinnlAC introduces a 
grid of magnitude and color shifts to the observed CMD. 
The values of the mean x 2 i n each node of the grid are 
analyzed to determine the best solution and the associ- 
ated shift from the adopted transformation to the abso- 
lute plane. In the process, the sets of 24 solutions were 
calculated 39 times, with different shifts, for a total of 
936 solutions. This computationally expensive task was 
perfor med using the Cond or workload management sys- 
ten0 (|Thain et al.1 120051 ) available at the Institute 'de 
Astrofisica de Canarias. An initial grid with 25 posi- 
tions was common to the two photometry sets. A second 
grid with 14 positions is used to refine the search for 
the best \ 2 . In the case of the DAOPHOT photometry, 
the minimum was confirmed at the position (S co i ; 5 mag ) 

= (0.0 ; 0.0), xl, min = 2 - 63 ' while in the case of the 
DOLPHOT photometry the minimum is located at (8 co i 
; S mag ) = (+0.06, 0.15), x%min = 2 - 57 - In the S31 we 
describe the details of the two solutions calculated at the 
xi, min position. 

4.2. The Match method 

The principles of the MATCH method of m easuring 
SFHs are presented in iDolphinl j2000a) and iDolphinl 
(2002). In this work, to appropriately compare the syn- 
thetic and observed CMDs, we converted them into Hess 
diagrams binning by 0.10 mags in Mf 475W and 0.05 mags 
in Mf^w — Mfsiaw- For consistency, we chose pa- 
rameters to create synthetic CMDs that closely matched 
those used in the IAC-pop code presented in this paper: 

17 1 http ://iac-star.iac.es\ 

18 We tested six differen t values, from / 3 =0 to 1 . A more detailed 
discussion can be found in Monclli et al. (2010T]). 

19 The value adopted in the SFH derivation, 24 .77, is slightly 
differ ent from the final value of 24.74 ± 0.12 given in DBernard et al.l 
12003) . However, such a small difference has negligible impact on 
the derived SFH. 

20 http:/ /www. cs.wisc.edu/condor/ 



• a power law IMF with a = 2.3 from 0.1 to 100 M ; 

• a binary fraction (3 — 0.40; 

• a distance modulus of 24.77, Af a7^w = 0.122; 

• the stellar evolution libraries of IMarigo et a l. (2008), 
which are the basic Padova/Girardi models with updated 
AGB evolutionary sequences. 

The synthetic CMDs were populated with stars over a 
time range of logi = 6.6 to logi = 10.20, with a uniform 
bin size of 0.1. Furthermore, the program was allowed 
to solve for the best fit metallicity per time bin, drawing 
from a range of [M/H] = —2.3 to 0.1, where M canoni- 
cally represents all metals. The depth of the photometry 
used for the SFH was equal to 50% completeness in both 
M_F475w and Mpsi4W- 

To quantify the accuracy of the resultant SFH, we 
tested for both statistical errors and modeling uncertain- 
ties. To simulate possible zero point discrepancies be- 
tween isochrones and data, we constructed SFHs adopt- 
ing small offsets in distance and extinction from the best 
fit values. The rms scatter between those solutions is 
a reasonable proxy for uncertainty in the stellar evolu- 
tion models and/or photometric zero-points. Statistical 
errors were accounted for by solving fifty random realiza- 
tions of the best fit SFH. Final error bars are calculated 
summing in quadrature the statistical and systematic er- 
rors. 

4.3. The COLE method 

The third method we applied to measure the SFH of 
Tucana is a simulated annealing algorithm applied to our 
DOLPHOT photometry and the Padova/Girardi library. 
The co de is the same a s that applied to the galaxies 
Leo A (ICole et al.l 120071) and IC 1613 (ISkillman et al.l 
|2003| ). and details of the method can be found in those 
references. 

In particular, the isochrone tables ([Marigo et al.l[2008l ) 
for each (age, metallicity) pair are shifted by the appro- 
priate distance and reddening values, and transformed 
into a discrete color-magnitude-density distribution g by 
convolution with: 

• the initial mass function from iChabrierl (|2003l ): 

• the results of the artificial star tests simulating obser- 
vational uncertainties and incompleteness; 

• the binary star frequency and mass ratio distribution 
from the solar neighborhood, namely 35% of stars are 
single, 46% are parametrized as "wide" binaries, i.e., the 
secondary mass is uncorrelated with the primary mass, 
and is drawn from the same IMF as the primary, and fi- 
nally 19% of stars are parameterized as "close" binaries, 
i.e., the probability distribution function of the mass ra- 
tio is flat. 

The history of the galaxy is then divided into dis- 
crete age bins with approximately constant logarithmic 
spacing. We used 10 time bins ranging from 8.00 < 
log(age/Gyr) < 10.13. Isochrones evenly spaced by 0.2 
dex from —2.3 < [M/H] < —0.9 were used, considering 
only scaled-solar abundances but with no further con- 
straints on the age-metallicity relation. The SFH is then 
modeled as the linear combination of ip(t,Z) that best 
matches the data by summing the values of g at each 
point in the Hess diagram. The CMD was discretized 
into bins of dimension 0.10x0.20 mag, and it is the den- 
sity distribution in this Hess diagram that was fit to the 
models. The error bars are determined by perturbing 
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Fig. 5. — 3-D representation showing the SFH of Tucana, derived 
with the DAOPHOT (top) and the DOLPHOT (bottom) photom- 
etry. Both show the average of 24 solutions calculated at the min- 
imum x 2 in the Smag-Scol grid. 

each component of the best-fit SFH in turn and re-solving 
for a new best fit with the perturbed component held 
fixed. 

4.4. Comparison of the different solutions 

We verified the consistency of the results obtained with 
the three methods, observing good agreement as shown 
in the previous papers. A comparison of the results 
obtained using different photometry/library/SFH-code 
combinations is shown in Figure^! which presents, in the 
panels a) and 6), the cumulative mass fraction as a func- 
tion of time for different combinations of photometry, li- 
brary, and method. In particular, the panel a)compares 
the solutions from the DAOPHOT and DOLPHOT pho- 
tometry derived using the BaSTI and Padova stellar evo- 
lution library. Panel 6) shows the solution derived with 
the three SFH codes, for the DOLPHOT photometry 
coupled with the Padova library. The two plots disclose a 
generally very good agreement between the different solu- 
tions. Panels c) and d) compare the i/>(t) calculated with 
the three SFH algorithms. It has to be stressed that the 



TABLE 1 

Integrated quantities derived for the Tucana dSph. 





DAOPHOT 


DOLPHOT 


fip(t)dt [1O 6 M ] 

< ip(t) > [10- 7 M Q j/r- 1 pc- 2 ] 

< age > [10 10 yr] 

< [Fe/H] > [10- 4 dex] 


(2.52±0.02) 
(1.25±0.01) 
(1.15±0.02) 
(5.78±0.83) 


(2.62±0.03) 
(1.30±0.01) 
(1.17±0.01) 
(5.63±0.82) 



three methods use different IMFs, and the comparison 
of the derived ip(t) should be done carefully. Therefore, 
panels c) and d) compare the three ^(t) before and after 
applying a scaling factor to the MATCH and COLE so- 
lutions, in order to take into account the different IMFs. 
Note the very good agreement in the main features: the 
position of the main peak at the oldest ages, and the sub- 
sequent decreasing trend. Both the MATCH and COLE 
methods give some small residual star formation between 
5 and 8 Gyr ago, which is not present in the IAC-pop so- 
lution, but still the three solutions are consistent within 
the error bars. 

4.5. The Tucana SFH 

Similarly to the p revious works (|Monelli et al.ll2010bl ; 
iHidalgo et al.l 12010ft , we discuss in detail the solution 
based on the IAC method applied to a model CMD cal- 
culated with the BaSTI library and, in particular, the 
average of the two solutions based on the two sets of 
photometry. 

We present two different representations of both so- 
lutions. Fig. [5] shows a 3D view of the best solutions, 
calculated as the run ning average of the 2 4 individual so- 
lutions at the xlmin (|Hidalgo et al.1 12010). This is a con- 
venient way to have an immediate, general visual ization 
of th e stellar content and evolution of a system (Hodge 
1989), since it includes the information on the star for- 
mation rate as a function of both age and metallicity, 
and the age and metallicity distributions as projections 
of the first. Fig. [5] shows that the bulk of star forma- 
tion occurred in Tucana at a very early time and with a 
limited metallicity range. Only a small percentage of the 
star formation inferred in the solution occurred at epochs 
younger than 10 Gyr ago and with metallicity Z > 0.001. 
Another feature that appears in both solutions is a small 
population of intermediate-age, low metallicity stars. We 
interpret this as a population of blue stragglers, because 
the inferred metallicity is in good agreement with that of 
the oldest population in Tucana, and it does not follow 
the general age-metallicity relation of the bulk of stars. 
T he same feature was pre sent in the Cetus SFH presented 
bv lMonelli et a l. (2010bJ), and was interpreted similarly. 
We refer to that reference for further discussion. 

A complementary view is shown in Fig. [5J where the 
star formation rate, the chemical enrichment law and 
the cumulative mass fraction are plotted as functions 
of time. The thick continuous and dashed lines show 
the DAOPHOT and DOLPHOT solutions, respectively, 
while the thin lines represent the corresponding error 
bars. These have been calculated following the prescrip- 
tion given in IHidalgo et al.l ((20101 ) . Note the good agree- 
ment between the solutions based on the two different 
sets of photometry. Both recover the highest peak as a 
flat feature at the oldest ages, t > 12.5 Gyr ago. The 




Fig. 6. — Summary of the Tucana SFH results. The three panels show, from top to bottom, tp(t), the age-metallicity relation and 
the cumulative mass fraction. We show the results with the BaSTI library for both the DAOPHOT (dashed line) and the DOLPHOT 
(continuous) photometry, and the average of the two (thick cyan line). The dashed part of the cyan line representing the average of the 
two sets of photometry in the middle panel indicates what we believe to be the presence of BSs rather than a low level of intermediate-age 
and young SF. The thin lines repre sent the error bars, w hile the vertical lines in the upper panel mark redshifts z = 15 and 6 — the epochs 
corresponding to the reionization. (Bouwcns ct al. 2007]). The dottcd-dashed lines in the central panel mark the epochs when 10%, 50% 
and 90% of the stellar mass was formed. 



steepness of the decline is marginally different, and both 
solutions give very similar ages at which the star for- 
mation completely stopped in Tucana, around 8-9 Gyr 
ago. There is also excellent agreement in the small peak 
detected at ~ 4 Gyr ago, related to the BS population. 
An excellent correspondence is also observed in the de- 
rived metallicity laws, represented in the middle panel. 
A steady increase of [M/H] is observed, from [M/H]«-1.6 
at very early times to [M/H]«— 1.1 at the end of the star 
formation activity. After that, the shape of Z(t) is nec- 
essarily uncertain given the small number of stars from 
which it is derived and, therefore, it is not shown in the 
plot. 

Since the small differences in ip(t) an d in the chemi- 
cal evolution law from the two photometry sets are well 
within the error bars and because we do not have any 
reason to prefer one photometry set over the other, it 
seems reasonable to adopt the average of the two solu- 
tions as the Tucana SFH, and the differences between the 
two as an indication of external errors. The thick cyan 
lines in Fig. |6] represent this average, which we adopt as 
the final solution for the Tucana SFH. 

The upper panel of Fig. [6] discloses that no star for- 
mation occurred in Tucana in the last » 8-9 Gyr. This 



statement can be further quantified by quoting the ages 
when we estimate that 10%, 50% and 95% of the stellar 
mass was formed (s=s 13.2, 12.1 and 9.7 Gyr ago; these are 
represented as vertical lines in the middle panel of Fig. 
[p. This result dem onstrates that Tucana, like Cetus 
(pvlonelh et al.ll2010bD . is an outlier of the morphology- 
density relation, because it is a purely old system located 
at a large distance from both the giant galaxies of the LG. 

Fig. presents a comparison between the observed 
DAOPHOT CMD and a best-fit CMD (left and central 
panels), represented as Hess diagrams. The bundles used 
to derive the solution have been superimposed on the ob- 
served CMD. We find satisfactory agreement in all the 
main evolutionary phases. As far as the differences in the 
HB morphology are concerned, it is worth noting that we 
have not tried to properly reproduce this observational 
feature. Due to the strong dependence of the HB mor- 
phology on the mass- loss processes along the RGB, it is 
clear that a better agreement between the observational 
data and the synthetic CMD could be obtained with a 
better knowledge of the average mass loss efficiency and 
its spread. The right panel shows the residual Hess dia- 
gram, in units of Poisson error, which supports the good 
agreement between the observed and the best fit CMD. 
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Fig. 7. — Comparison of the observed (left) and best fit (center) CMDs. The latter was derived by extracting random stars from the 
model CMD used to derive the solution, in such a way that each simple population contributes proportionally to the calculated SFR. The 
black lines in the left panel mark the position of the four bundles adopted to derive the solution. The right panel shows the residual Hess 
diagram (calculated as observed— model) , shown in units of Poisson errors. The Figure discloses an overall general agreement between the 
observed and the best-fit CMDs, particularly in the whole TO region used to calculate the SFH. The discrepancy at the faintest magnitude 
is due to the limiting magnitude of the model CMD. Note that such a discrepancy does not affect the portion of the CMD included in the 
bundles. The differences in the HB morphology are expected, since small differences in the adopted mass-loss prescription along the RGB 
have a big impact on the HB morphology. 



The mean quantities derived (star formation rate, age, 
and metallicity) , together with the integral of the ^(t), 
are summarized in Table [TJ 

4.6. The short first episode of star formation in 
Tucana: are we seeing the effects of reionization? 

The description presented in the previous section cor- 
responds to the SFH of Tucana as derived using the algo- 
rithms discussed above. This solution is the convolution 
of the actual Tucana SFH with some smearing produced 
by the combination of observational errors and the limi- 
tatio ns of the SFH recovery method (|Aparicio fc Hidalgo! 
12001 . These are related in part to intrinsic limitations 
of the method, such as the small differences in position 
in the MSTO of old stars of ages within 2-3 Gyr. In the 
following, we will try to constrain further the features of 
the actual, underlying Tucana SFH through some tests 
with mock stellar populations. 

In particular, to address the question: "Did cosmic 
reionization affect the SFH of Tucana?" we designed a 
specific test with mock populations. It is generally as- 
sumed that the Universe was completely reionized by 
z w 6, equivalent to a look-back time of ~ 12.8 Gyr 
(Bouwens et al 2007). Therefore, we created three mock 
stellar populations with constant ip(t) in the age ranges 



13.3 < t < 12.8 Gyr, 12.8 < t < 12.3 Gyr, and 12.3 
< 11.8 Gyr agdZJ We also assumed a mean metallicity 
and a metallicity dispersion similar to what is inferred 
for Tucana. We simulated the observational errors of the 
DAOPHOT photometry, and recovered the SFH of the 
mock galaxies with the same prescription as for the real 
galaxy. In order to minimize the fluctuations associated 
with the creation of the synthetic CMDs, we repeated 
each test five times. Figure [51 left panel, shows the mean 
of the five resulting ip(t) for each mock population, com- 
pared with that of Tucana, shown as a thick black line 
with the corresponding error bars. Note that the fact of 
repeating the experiment five times reduces the impact 
of the Poisson errors, reflected in the smaller error bars 
of the recovered mock ip(t). The time extension of the 
three input bursts is shown schematically at the bottom 
of the figure. The three recovered functions ip(t) show 
similar general behavior, but the main peak is recovered 
younger and lower, and its width larger, for decreasing 
age of the input ip(t). The comparison with the Tucana 
SFH leads to two main conclusions: 
• despite the fact that the best agreement at the oldest 
epoch is obtained with the intermediate-age mock pop- 
ulation, the case of the oldest mock population is also 
compatible, within the error bars, with the SFH observed 
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Fig. 8. — Results of the tests performed recovering the SFH of two mock populations. The left panel compares the Tucana ^((t) with 
that of three mock bursts at different ages. The central and right panels compare it with Gaussian profile ip(t) laws of different a and mean 
age. These comparisons disclose that the bulk of the star formation in Tucana occurred between 13.5 and 12 Gyr ago, suggesting that we 
can not put firm constraints on the effect of the reionization in this galaxy. Nevertheless, these tests strongly suggest that Tucana could 
not undergo strong star formation at epochs younger than 12 Gyr ago, although some residual star formation is still compatible with the 
observed $(t). 



for Tucana. This suggests that, due to the limited age 
resolution, we can not draw definitive conclusions con- 
cerning the effect of the reionization on Tucana. In fact, 
it is equally acceptable that Tucana started forming stars 
either before or after the Universe was completely reion- 
ized, more than 12.8 Gyr ago. Moreover, note that if the 
star formation started before 12.8 Gyr ago, the derived 
solution may indicate that the star formation extends 
somewhat toward younger epochs. 

• The youngest mock population is only marginally con- 
sistent with the Tucana solution at all ages, strongly sug- 
gesting the the bulk of star formation in this galaxy could 
not occur at epochs younger than ~12 Gyr ago. How- 
ever, it is still possible that some residual star formation 
was active at younger ages. 

21 Note that the age derived through the SFH recovery depend 
on the assumed set of theoretical models. However, it has been 
shown that the libraries adopted in this work give similar ages 



Note that the mock populations presented in Figure 
[8] are characterized by a relatively short episode of star 
formation, lasting 0.5 Gyr. Nevertheless, the recovered 
ip(t) shows a significantly longer duration, comparable to 
that of Tucana. This suggests that the actual duration of 
the main star formation event in Tucana was significantly 
shorter than the calculated one. 

To further investigate this point, we simulated more 
mock populations, with Gaussian input ?/>(t) peaked at 
13 Gyr ago and with age dispersion a = 0.1, 1.0 and 
1.5 Gyr. The first two have also been repeated assum- 
ing an age peak at 12 Gyr ago. The Gaussian tp(t) was 
truncated at the oldest age we are considering, 13.5 Gyr, 
i.e., no stars older than 13.5 Gyr were allowed in the 



IjGallar t ct al. 2005), and also that the ages derived for the Galactic 
globular clusters are in excellent agre ement with the age of th e 
Universe from the WMAP experiment ( Marin- Franch ct al. 2009). 
This ensures a meaningful comparison between the solution SFH 
in the same cosmological framework. 
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14 13 12 11 10 9 8 7 
Look-back time (Gyr) 

Fig. 9. — Comparison of the ip(t) {top), age-metallicity relations 
(middle), and cumulative mass functions of Tucana and Cetus. The 
epoch of peak star formation in Tucana is ~ 1.0 Gyr earlier than 
in Cetus. 



mock populations. The results are shown in Figure [5J 
central and right panels. As in the previous test, we re- 
peated each case hve times and show the mean recovered 
ip(t). The comparison in the central panel shows that the 
observed peak at the oldest ages is best reproduced by 
the (7=1.0 Gyr Gaussian profile, but it is still compatible 
with the narrowest one. 

On the other hand, the narrowest Gaussian profile does 
not properly represent the shape of the Tucana ip(t), sug- 
gesting that the first event was somewhat longer than 
this test case. Similarly, the peak of the widest Gaus- 
sian profile appears too low and flat, and is not com- 
patible with the observed one. The comparison with the 
Gaussian ip(t) peaked at 12 Gyr ago, shown in the right 
panel, indicates that the bulk of the star formation of 
these mock populations is too young and is not compat- 
ible with the Tucana ip(t). Note that, in equivalent te sts 
performed for the Cetus dSph in Monelli et al. ( 2010bl). a 
mock Gaussian population peaked at this age was a good 
representation of the observed ip(t). This suggests that, 
despite the apparent similarities of the calculated SFH 
of these two galaxies, we are able to disentangle subtle 
age differences between the two (see [J5]). 

We stress that the case of Tucana is possibly the most 
difficult among the LCID galaxies. The combination 
of both intrinsic effects (it is the most distant galaxy 
in our sample, and the oldest one) and observational 
ones (the shortest exposure times due to the South At- 
lantic Anomaly) , contribute to increase the complexity of 
the problem and the loss of age resolution at the oldest 
epochs. 



4.7. Summary 

The SFH of Tucana can be summarized as follows. Tu- 
cana experienced a very early and short event of star for- 
mation, peaked at ages t > 12.5 Gyr ago. We measure 
a steady decrease in the star formation rate until star 
formation completely stopped ~ 9 Gyr ago. We showed 
that the measured duration is overestimated and that 
the actual event is compatible with a Gaussian ip(t) with 
a between < 1.0 Gyr. The apparent residual star forma- 
tion at intermediate ages is interpreted as a population 
of blue stragglers. The tests performed with mock popu- 
lations suggest that the epoch of the main episode of star 
formation is compatible with the end of the reionization 
era, and that it can not have occurred later than 12 Gyr 
ago. This is different from what found for Cetus, where 
the majority of stars formed at least 1 Gyr later. We 
will present an in-depth comparison of the two galaxies 
in the next section. 

5. TUCANA AND CETUS 

Tucana and Cetus share the noticeable characteristic 
of being among the most isolated dSph galaxies in the 
LG. Cetus is located at more than 650 Kpc from both the 
Milky Way and M31. Tucana is even more distant, at 900 
Kpc from the Milky Way and 1.3 Mpc from M31. The 
actual position and orbit are crucial to understand inter- 
actions that may have occurred in the past. Their reced- 
ing velocities with respect of the barycen ter of the LG 
(jFraternali et all l2009t iLewis et all 12007ft suggest that 
both galaxies may already have experienced at least one 
passage through the innermost regions of the LG and, in 
the case of Cetus, that it is close to the apocenter of its 
orbit. 

Therefore, a detailed comparison of their properties 
can be useful to understand the origin of the differences 
and similarities in their SFHs. 

5.1. Physical parameters 

Tab. [5] gives a summary of the intrinsic properties 
of Cetus and Tucana. Cetus is the larger galaxy, both 
in terms of baryonic mass and physical size. Tucana 
presents the mass (1.81 x 10 6 M e baryonic mass, Hidalgo 
et al. in prep) and size (exponential scale length a ~ 120 
pc) typical of the smallest dSphs in the LG, similar to 
Sculptor, Leo II and Ursa Minor. On the other hand, Ce- 
tus is similar to the larger nearby dSphs, with mass and 
size comparable to Carina or And II, with a total bary- 
onic mass of 3.94 x 10 6 Mq and a ~ 220 pc. Despite their 
large distances from both the Milky Way and M31, both 
Tucana and Cetus are similar to the oldest systems like 
Draco, Ursa Minor or Sculptor. However, the velocity 
dispersion estimated for the two galaxies (a ~ 17 'kms~ l ) 
is significantly higher than the measured dispersion for 
the nearby dSphs (< lOfcms -1 ). Assuming that the cir- 
cular velocity of the halo v c is v c = y/[3)a, and the scal- 
ing suggested for su b-halos by cosmological simulations 
(see iKravtsovl 120101 ) and confirmed also by simulations 
of tidally perturbed satelli tes with dark matter and stars 
(|Klimentowski et alJl2007j h one obtains v c ~ 29 km/s for 
the halo of Tucana, which would correspond to a virial 
mass of ~ 4 x 1Q 9 M P ,, higher than that of dSph s atel- 
lites of the MW rtLokasl [200l IStrieari et alll2010ft but 
comparable to the mass estimated for their progenitors 
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Fig. 10. — Comparison of the Tucana (left) and Cetus (right) i/>(t) relations with those of Gaussian mock populations giving comparable 
solutions. The inferred underlying SFHs of Tucana and Cetus support the idea of an earlier and stronger first episode in Tucana. 



be fore they were accreted by the Milky Way or M 31 ha- 
los (|Maverll20ia see the reviews bv lKravtsovll20irj, and). 
Since Tucana and Cetus have been only weakly affected 
by tides the estimated mass should indeed be close to 
their mass before infall into the primary galaxies. With 
such a high halo mass the resulting M/L ratio would 
be several hundred in solar units. However, both show 
hints of rotation, and taking this into account would re- 
duce these numbers. Moreover, the faintness of the ob- 
served stars, together with the low S/N of the available 
low-resolution Ca-triplet spectra might also have played 
a role in overestimating the velocity dispersion. It is 
nonetheless conceivable that M/L is at least as high as 
in the most dark-matter-dominated among the classical 
dSphs, such as Draco, Ursa Minor and And II. 

5.2. The oldest stars as a benchmark for the 
reionization effects 

Fig. [S] compares the derived properties for Tucana 
(black line) and Cetus (grey line). The upper panel shows 
the calculated ip(t) as the average from the two pho- 
tometry sets. All the experiments performed to test the 
IAC-pop code (lAparicio fc Hidalgol 120091 : iHidalgo etall 
l2010t lMonelli et all 2010bO . as well as the tests with mock 
galaxies presented in £)4.6[ support the reliability of the 



solutions, in relation to the determination of the age of 
the peak star format ion (the typical prec ision is of the 
order of 0.5 Gyr, see lMonelli et al.ll2010bT) . including at 
epochs > 10 Gyr. Moreover, we stress that we are com- 
paring two SFHs obtained using the same stellar evo- 
lution library, and the average resulting from the two 
photometry sets. This minimizes the possible systematic 
effects, giving more confidence to the estimate of the rel- 
ative differences between the main SFH features of the 
two galaxies. 

The analysis of the ^(t) shows that the evolution of 
the two galaxies was different. Tucana experienced its 
strongest level of star formation at the earliest epochs, 
with the initial peak followed by a steady decrease in the 
star formatio n rate. On the other hand, the star forma- 
tion in Cetus (jMonelli et alJl201QbT ) might have started at 
the same epoch, but with the rate increasing until reach- 
ing the main peak ~ 1.0 Gyr later. Therefore, there is a 
strong indication that Cetus formed the majority of its 
stellar mass significantly after Tucana, as summarized in 
Tab. [2] with additional information such as £50 and the 
mean age. 

Note that we have observed most of Tucana's body, 
while in the case of Cetus, the very central region of 
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the galaxy (within K gc < 0.5r c ) was not contained 
in our ACS data. This may affect somewhat the re- 
sulting SFH, in view of the stellar population gradi- 
ents that are present in these galaxies (which will be 
discussed in detail in Hidalgo et al., in prep.). How- 
ever, in this study we found — in agreement with pre- 
vious works ([Gallagher et al.lll99§ iHarbeck et al.ll2001l; 
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20031: IB attagli a et al.l 



20061: IBernard et alj 



20081: iLee et al.l 120091 ) — that the stellar population gra 



dients are in the sense that the population is on average 
younger in the central part. This implies that the former 
picture of the differences between the SFH of Tucana and 
Cetus would unlikely change (or if changing, the differ- 
ences would only increase) if we could compare with the 
global SFH of Cetus. 

The central panel of Figure [9] shows the comparison of 
the two age-metallicity relations. The thin lines at ei- 
ther side of each relation indicate the 1-a dispersion of 
the metallicity distribution at each age; the whole metal- 
licity range at each age is much larger. The run of mean 
metallicity with age is similar in both galaxies, with the 
mean metallicity steadily increasing as a function of time. 
However, Tucana appears more metal-rich than Cetus, 
at fixed age, in contrast w ith the metallicity- luminosity 
relation of dSph galaxies (jGrebel et al.l 12003). but still 
consistent with the scatter of that relation. Note that 
the higher mean metallicity in Tucana does not imply a 
difference in the metallicity of the most metal-poor pop- 
ulations, which is the same in both galaxies and consis- 
tent with Z=0.0001. The higher mean metallicity derived 
for Tucana may also be connected to the stronger initial 
event of star formation, which enriched the medium in a 
short time scale that we are not able to resolve. 

The lower panel of Figure [9] compares the cumulative 
mass fraction as a function of time. Tucana presents 
a steeper curve, reflecting the discussed stronger initial 
activity. 

As discuss e d in iHidalgo et all ()201 0) and in 
iMonelli et al] (|2010bl) . the effect on the derived 
SFH of our limited age resolution at old ages, is mainly 
to broaden the main features, and to artificially extend 
the limits of the epoch(s) of star formation. The subtle 
differences between the Tucana and Cetus SFHs may be 
best visualized in terms of the possible true, underlying 
SFHs that we have discussed for both galaxies, using 
tests with mock stellar populations. Figure [10] summa- 
rizes the comparison of the tp(t) retrieved for Tucana 
and Cetus with those of the compatible Gaussian mock 
populations. Despite being only qualitative, this plot 
clearly discloses that the star formation in Tucana is 
expected to have been significantly stronger at the 
oldest ages than in the case of Cetus and that, since 
the narrowest events seem to be compatible with the 
SFH of both galaxies, even the main star formation 
event may have occurred at almost non-overlapping 
epochs. Despite the uncertainties possibly affecting the 
absolute age scale, and the ambiguity concerning the 
duration of the main event in Tucana highlighted in 
§4.6| the relative difference of 1 Gyr calculated between 
the inferred peak of the V'(t) of Tucana and that of 
Cetus strongly reinfo rces the conclusions presented in 
IMonelli et al.1 (|2010b| ). where we argued that the reion- 
ization had limited, if any, effects on the star formation 
of Cetus. This is expected if the halo masses of these 



dwar fs are > 1O 8 M (jSusa fc Umemural 120041: iKravtsovl 
l2OT0l) . which can be easily argued even if their velocity 
dispersion is somewhat overestimated. We can rule out 
the possibility that Cetus experienced the majority of 
its star formation at epochs earlier than the end of the 
reionization. This conclusion is solid and is further 
sustained by the consistent results obtained using 
different photometry sets, stellar evolution libraries and 
SFH codes adopted for the LCID project. 

5.3. Morphology of the HB and RR Lyrae stars 
properties 

Is it possible to link the resulting SFHs with other 
independent evidence? A comparison of the Tucana and 
Cetus CMDs is shown in Fig. [2j At a first glance, the 
striking similarities seen in the shape of the TO region 
and in the slope and color width of the RGB suggest a 
similar age range and spread of metal content. 

An obvious difference between the two shows up in 
the morphology of the HB. The HB in Tucana presents a 
complex morphology, with both the blue and the red side 
well populated, the red side also presenting a split into 
two paral lel sequences. Cetus presents a redder HB mor - 
phology (jSaraiedini et al.1 I2002t IMonelli et all I2010bl) . 
but nonetheless the blue side is populated, even if by 
a many fewer stars. The color of the bluest Cetus HB 
stars is the same as for Tucana. 

Therefore, the color limits of the HB are similar in the 
two galaxies, at both extremes. If we assume that, in 
complex systems like these, the prolonged star formation 
and the associated chemical evolution drive such a color 
spread, then we can expect that the blue side is popu- 
lated by the oldest, most metal-poor stars. Therefore, 
the higher number of stars in this region of the Tucana 
CMD suggests a much stronger star formation event at 
older epochs and/or lower metallicities than for Cetus. 
However, the presence of a few blue HB stars in the lat- 
ter supports the idea that some amount of very old pop- 
ulation s present in both. These considerations are in 
excellent agreement with the derived SFHs (see Figures 

The same scenario is also nic ely supported by the 
properties of the RR Lyrae stars. IBernard et al.l (2008) 
showed the presence of two sub-populations of RR Lyrae 
stars in Tucana, characterized by different mean lumi- 
nosities, different mean periods, different slopes of the 
period-amplitude relation, and different radial distribu- 
tions. The interpretation of this evidence is that the two 
sub-populations have different properties in terms of age 
and/or metallicity. Note that the fainter sample, pre- 
sumably the younger and more metal-rich, is also the 
more centrally concentrated. A similar dichotomy was 
not detected in Cetus, where the sample of RR Lyrae 
presents homogeneous properties. 

The presence of two sub-populations of old stars in Tu- 
cana, supported by both the split in the red HB and the 
properties of the RR Lyrae stars, implies that that the 
first generation is numerous enough to show up clearly 
and that Tucana was able to produce a second generation 
of stars with distinct chemical properties due to the self- 
enrichment from the previous event. The time resolution 
achievable with the present data is not enough to distin- 
guish an episodic from a continuous ?/>(t). Nevertheless, 
the chemical enrichment seems to have occurred on very 
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TABLE 2 

Estimates of the basic parameters for the two dSph galaxies Tucana and Cetus 





Ti ipQTin 




D [Kpc] 


887 ± 50 d 


780 ± 43 d 


(m-M)o 


24.74 ±0.12 d 


24.46 ±0.12 d 


My [mag] 


-9.55 ±0.27° 


-10.1 a ;-11.3 b 


Lv [L0] 


5.2 x 10 5e 


8.7 x 10 5a ;2.6 x 10 6b 


M/L [M0/L0] 


105±|| c 


~ 70 f 


Total Mass [M0] 


~ 6 x 10 7c 


~ 6 x 10 7f 


Stellar Mass [M©] 


1.81 x 10 6 


3.94 x 10 6 


a['(pc)] 
Rc ['(pc)] 


0.48 ± 0.08(124) c ; 0.51(132)8 


1.0 ± XX(230) g ; 1.59 ± 0.05 b 


0.59(152) d ; 0.66 ± 0.1(170 ± 25) c 


1.5 ± 0.1(340) a ; 1.3 ± 0.1(295) b 


Rt ['(pc)] 


3.45(890) c 


4.8 ± 0.2(1089 ± 45) a ; 32 ± 6.5(7261 ± 1475) b 


c 


0.72 c 


0.5 a ;1.4 b 


Vo,v (obs)[mag] 


25.05 ±0.06° 


25.0 b 


Vhel [km/s] 


+ 194.0 ±4.3 C 


-87 ± 2 f 


vasR [km/s] 


+98.9 C 


-25 f 


VLGSR [km/s] 


+73.3 C 


+14 f 


a [km/s] 


1 7 A -3.1c 


17±2(14H-8) f 


vrot [km/s] 


~ 16 c 


7.7 ± 1.2 f 


<age> [Gyr] 


11.6±0.1 h 


11.3 ± 0.2S 


Tio [Gyr] 


13. 2 h 


13. I s 


T50 [Gyr] 


12. l h 


11.68 


T95 [Gyr] 


9.7 h 


9.6 s 



a IWhiting et al.H199Sl; 

b McConnachic k. Irwin 2006; 

c Fratornali ct al. 2009; 

d Bernard ct al. 2009; 

c Savianc ct al. 1996; 

f ILewis et a.1.1120071; 

g Monclli ct al. 2010b; 

b this work. 



short time scale. The fact that the mean metallicity is 
relatively high even at the oldest epochs (see Figure 
central panel), indicates that the enrichment occurred 
very early. Therefore, the derived SFH, dominated by a 
single event, is still compatible with two populations of 
old stars identified through other evidence. 

5.4. The RGB bump 

The RGB bump is a well-define d feature in the CMP 
of both galaxies (see Figure land iMonelli et afll2010af ). 
In particular, Tucana presents two bumps, a bluer and 
brighter one, well separated from a second which is 
redder and fainter. Cetus only shows one well popu- 
lated bump, l ocated on the red side of the RGB. In 
IMonelli et all (|2010aD . we discussed how these observed 
features can be interpreted using the SFHs of the two 
galaxies, and we put constraints on the age and the 
metallicity of the populations which produce them. More 
metal-poor stars (Z<0.0005) produce the brighter bump 
in Tucana, while the fainter, which is also more popu- 
lated, consists of slightly more metal-rich stars (0.0005 < 
Z < 0.001). In Cetus, the dominant population (0.0003 < 
Z < 0.0006) produces the single bump. The RGB bumps, 
therefore, also reflect the previously discussed differences 
in the Tucana and Cetus SFHs, with an older, more metal 
poor population in Tucana numerous enough to show in 
a separate RGB bump in the blue part of the RGB. 

5.5. Summary 



All these pieces of information give a coherent picture 
of the evolution of Cetus and Tucana and of the origin of 
the subtle differences we have detected. The SFH of Tu- 
cana is compatible with a very old and brief event. The 
presence of two sub-populations of old stars in Tucana, 
supported both by the split in the red HB, the proper- 
ties of the RR Lyrae stars and the double RGB-bump, 
implies that Tucana was able to produce a second gen- 
eration of stars with distinct chemical properties due to 
sclf-cnrichmcnt from the first stars. The important frac- 
tion of very metal-poor stars in Tucana indicates a high 
efficiency of star formation at the earliest stages. We 
cannot distinguish, with the present data, between two 
discrete events or a continuous process. 

On the other hand, Cetus started forming stars most 
likely at the same epoch, as indicated by the bluest HB 
stars. However, Cetus seems to have had a more gradual 
evolution, reaching its peak star formation 1 Gyr after 
Tucana, and producing a dominant population with more 
homogeneous metallicity. It contains a minority of old, 
metal poor stars, like the oldest stars in Tucana, and a 
majority of slightly younger stars (and possibly some- 
what metal enriched with respect to the oldest). The 
different mean metallicities of Cetus and Tucana can be 
possibly explained by the intense initial episode of star 
formation that Tucana experienced at older epochs, that 
enriched the medium causing the higher metallicity of 
the subsequent generation of stars. 

According to the most widely accepted cosmological 
framework, the age of the Universe is ~13.7 Gyr, and the 
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best estimate of the re ionization epoch is at z ~ 10.9±1.4 
(|Komatsu et al.l [20091 ). Therefore, the most important 
conclusion we can draw from the comparison of the SFH 
of Cetus and Tucana is that the vast majority of the stars 
in Cetus were born after the end of the reionization, cur- 
rently assumed at z ~ 6. This confirms that the reioniza- 
tion can not have been the domina nt mechanism causin g 
the end of star formation in Cetus (|Monelh et al. 2010b). 
On the other hand, the timescale of the star formation in 
Tucana suggests that the reionization may have played 
an important role in this galaxy. 

6. DISCUSSION 

The SFH of the Tucana isolated dSph galaxy presented 
in this paper shows that the star formation was very ac- 
tive at very early epochs (> 12.5 Gyr ago) and over a 
short time scale (compatible with a Gaussian ip(t) with 
a in the range t < 1.0 Gyr). The tests we performed 
with mock galaxies show that the present data do not al- 
low enough time resolution to firmly correlate the epoch 
of the reionization with the end of the star formation in 
Tucana. The comparison with Cetus discloses that the 
latter undoubtedly experienced delayed activity, as com- 
pared with Tucana, and the vast majority of its stars 
were formed significantly later than the end of the reion- 
ization era. This would be compatible with the higher 
mass of Cetus, and with the f act that lower m ass systems 
are expected to form earlier (jKravtso v 2010). However, 
it is particularly intriguing that LGS3 is found to be less 
massive than Tucana, but nevertheless it formed 80% of 
its stars well after the end of the reionizatio n, with a peak 
star f ormation at the same age as in Cetus (iHidalgo et al.l 
f2010h . 

Given this, if we assume that the end of the star forma- 
tion was not predominantly induced by the reionization, 
what was the cause? In the following, we explore possi- 
ble alternative mechanisms, both internal and external, 
that might have affected the evolution of Tucana. 

6.1. The effect of SN explosions 

We will explore here the possibility that the star for- 
mation in Tucana was stopped by the removal of the gas 
by purely internal processes, such as the explosion of su- 
pernovae (SN e). Given the adopted IMF (assumed from 
Kroupa 2001), coupled with the accurate estimate of the 
total mass transformed into stars in the galaxy (obtained 
as the integral of the star formation in the main event, re- 
ported in Tabled]), it is possible to calculate the expected 
total number of type II SNe originating from the explo- 
sion of massive stars. At the metallicity typical of the 
Tucana population, we can assume that SNe II originate 
from stars with masses larger than f» 6.5M Q . Integrat- 
ing the IMF, and assuming 3.2 x 10 6 M Q to be the mass 
of stars formed as derived from the SFH, we estimate a 
total number of ^54,300 SNe II events. This total num- 
ber does not depend on the duration of the burst, but 
the number of SNs per unit time increases for decreasing 
duration of the star-formation event. For the two cases 
examined in tj4.6l we can assume a duration equal to the 
FWHM of the input Gaussian: 0.235 Gyr and 1.67 Gyr 
for the a —0.1 and 1.0 Gyr cases. Note that, in the latter 
case, the assumed duration is smaller than the theoretical 
value of 2.35 Gyr because we consider only stars younger 



than 13.5 Gyr. Assuming that each SN releases 10 51 erg, 
this implies a total luminosity of 7.33 x 10 39 ergs -1 and 
1.03 x 10 39 ergs -1 , for the 0.1 and 1.0 Gyr events, respec- 
tiveljO- If we compare these numbers with the model 
proposed by iMac Low fc Ferraral (|1999l see their Figure 
1), we can see that in both cases Tucana is in the so-called 
blowout regime, meaning that part of the gas must be lost 
due to the wind produced by SNe. This gives a general 
indication that gas loss, and especially the loss of metals 
directly connected with the SN wind s, are efficient mech- 
anisms in small systems like Tucana (IMac Low fc Ferraral 
flflOllFerrara fc Tolstovll2?M ISawala et al.H2010D .~ 

However, we have shown strong evidence that two stel- 
lar populations with different properties formed in Tu- 
cana on a short time scale, and that the younger popula- 
tion was enriched by the yield of the previous one. This 
has two important implications: i) Tucana was able to 
retain a sizable fraction of its gas after a first episode, 
and ii) the time between the two episodes of star for- 
mation was long enough to enrich and mix the medium. 
Given these two results, it is possible to estimate both 
the minimum time scale required to start enriching the 
medium, and amount of yield that has been produced 
by the first generation and recycled by the second one, 
and thus to estimate the amount of metals lost in an 
independent way. 

In the case of SNe II, the typical time scale is extremely 
short, of the order of few Myr. In the case of SNe la, this 
is related to the evolution of intermediate-mass stars in 
binary systems, which ends up with a degenerate CO core 
of mass insufficient for quiet ignition of the C burning. 
The mass limit of the star that produces a degenerate CO 
core has been shown to depend on the metallicity of the 
star (|Cassisi fc Castcllani 1993). For the assumed peak 
metallicity of the oldest population, the star with largest 
initial mass capable of producing a degenerate CO core 
is expected to be « 6.5M© star. The evolutionary time 
of such stars fixes a minimum time scale for the creation 
of SNe la of ~ 60 Myr. This suggests that, even in the 
scenario of the shortest Gaussian profile burst, SNe la 
were capable of contributing to the enrichment of the 
interstellar medium of Tucana. 

Second, it is possible to estimate the number of SNe re- 
quired to produce the iron that enriched the more metal- 
rich population in Tucana. The morphology of the HB 
suggests that the two populations of old stars are com- 
patible with mean metallicities Z = 0.0002±0.0001 and 
Z = 0.0006±0.0002, according to the comparison with 
theoretical ZAHBs and the SFH results. Based on the 
star counts on the HB, where the two populations can be 
more easily split, we find that each of the two episodes 
formed ~50% of the stars. Assuming the above metal- 
licities, and considering the total mass of stars formed, 
we can calculate the total mass of the metals contained 
in the younger population. In particular, assuming that 
each episode formed 1.6 xl0 6 Af Q , it turns out that 640 
Mq of heavy elements are recycled by the second genera- 

2 Note that this can be considered a conservative upper limit, 
since the total number of expected SN critically depends on the 
minimum stellar mass assumed to produce a SN event. This, in 
turn, strongly depends on the overshooting and the mass loss expe- 
rienced by the stars during their evolution: the effect of mass loss 
is to increase the mass of stars ending as a SN, therefore decreasing 
the number of events. 
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tion. Assuming the Solar mixture from lGrevesse &: Noelsl 
(1993), the amount of iron expected is ~ 46M . If we 
assume that only SNe II are involved, using the iron mass 
production rate from lNomoto et all (jl997b[ ). we estimate 
that ~ 300 SNe II are enough to produce this amount of 
iron. This estimate was derived assuming the iron mass 
production of a 13 Mq star. This value might change by 
a factor of 2-3 with different assumptions. We can now 
compare this number with the total expected number of 
SNe II derived before. If we take into account that stars 
more massive than ~ 25M Q do not enrich the medium, 
because they end up as neutron stars or black holes, we 
derive that, from a total of ^23,540 SNe II releasing 
iron during the first star formation event, only «1.3% 
of the iron mass is recycled, while the vast majority is 
lost, in agreement with the model bv lMac Low fc Ferraral 
(1999). On the other ext reme, since SNe la a re more 
efficient producers of iron (jNomoto et al.|[l997al) . we es- 
timate that only ~70 SNe la are required to produce the 
iron present in the younger population and that therefore 
an even larger fraction of gas is lost. 

6.2. The effect of the environment 

The SFH presented in this paper proves that Tucana, 
like Cetus, is an outlier in the morphology-density rela- 
tion. Despite its large distance from both the Milky Way 
and M31, Tucana is a purely old galaxy, similar to the 
oldest satellites of the Milky Way like Draco, Ursa Mi- 
nor, or Sculptor, with which it seems to share also a very 
high M/L. This poses interesting questions on the effec- 
tiveness of the environment in shaping the SFH. Why is 
a galaxy that spent most of its life in isolation so similar 
to galaxies of similar mass which are close satellites? The 
idea that close encounters with giant galaxies can influ- 
ence the evolution of d warf galaxies is a well-established 
prediction of mo dels (jMaver et all 120011 : iMayerl 120101 : 
lLokas et al.ll2010T) . However, at least 2-3 orbits on rela- 
tively small pericenters (< 50 kpc) are typically required 
to completely transform an initially gas-rich disky-like 
dwarf into a dSph, leaving a nearly isotropic stellar com- 
ponent suppo rted by stellar velocity dispersion ra ther 
than rotation (|Klimentowski et al.l[200l lLokasll2009f ) . A 
single passage on a relatively wide orbit (for an orbit with 
a 1:10 ratio between apocenter and pericenter, which 
is slightly larger than the mean but still common for 
sub-halos in cosmolo gical LCDM simulations— see, e.g., 
iDiemand et all 120071 and iKlimentowski et al.1 120101 the 
pericenter distance of Tucana must have been at R > 50 
kpc) would produce a remnant retain ing significant ro- 
tatio n and deviations from sphericity (Kazantzidi s et all 
l2010f) . While the residual rotation in Tucana argues in 
favor of the latter possibility, its structure and kinemat- 
ics are not firm enough yet to allow drawing conclusions 
of this sort. 

However, based on actual knowledge of the orbit of 
Tucana, it is concluded that it may have experienced at 
most one close encounter with a LG giant galaxy. In 
fact, the velocity of Tucana with respect to the barycen- 
ter o f the LG, vlgsr = +73.3fcms _1 (jFraternali et al.l 
2009) , suggests that a close encounter in the densest LG 
regions was possible, and in particular the velocity with 
respect to the Milky Way, vqsr = +98.9kms~ 1 , is com- 
patibl e with an encounter t hat occurred at t enc <10 Gyr 
ago ()Fraternali et al.l [2009) . This is a somewhat later 



than the main star-formation epoch identified in this 
paper, which appears to have ended between 12.5 and 
11.5 Gyr ago, but two events could have been nearly si- 
multaneous, given the large uncertainties. A single pas- 
sage at a distance of tens of kpc could be enough to 
strip most of the interstellar m edium of the galaxy by 
ram pressure (|Maver et al.ll2007h ■ provided that the coro- 
nal gas density around the proto-Milky Way was as high 
as suggested by observational constraints on its present 
state (~ 10~ atoms /cm 3 ). Heating by supernova feed- 
back and by the cosmic ionizing background, which re- 
mained high down to z ~ 2, would have rendered the gas 
more diffuse and pressure-supported, facilitating strip- 
ping by the combined action of ram pressure and tides 
(|Maver et all 120061: lMaverll2010l ). However, it should be 
taken into account that at such early epochs the size of 
the Milky Way and M31 was significantly smaller, so the 
effects of a close encounter might not have been as sig- 
nificant. Also, the structure of the coronal gas at such 
high redshift, when gal axies are undergo i ng rapid gas 
accre tion via cold flows (|Dekel et al.l 120091 : lAgertz et al.l 
2009), is unclear; while in principle gas densities around 
the primary could be even higher than at present as the 
filamentary flows are made of cold relatively dense gas, 
the halo might not be filled with gas everywhere as in a 
smooth corona. 

Moreover, if it turns out that Tucana has little 
or no residual rotation a complete transformation via 
tidal perturbations after only one pericenter passage 
would require quite specific conditions, based on grav- 
itational resonances o ccurring during the interaction 
(jD'Onghia et al.ll2009D . Another possibility is that Tu- 
cana was pushed to its high ly eccentric orbit as a result 
of a three-body interaction (S ales et al.1 12007). and thus 
could have had more than one close encounter with one 
of the primary galaxies. Finally, in cosmological simula- 
tions a few sub-halos undergo interactions and mergers 
with ot her sub-halos before being accreted by the pr i- 
maries (|Kravtsov et al.l 120041: IKlimentowski et al.ll2010D . 
which could also produce a dSph out of two interacting 
small disky dwarfs well outside the virial radius of the 
primary. 

On the other hand, the evidence that the properties of 
nearby and isolated dSphs present important similarities 
might indicate that interactions are not the dominant 
mechanism shaping these galaxies. For example, the su- 
pernova explosions from the first generation of stars con- 
tributed to the turbulence of the gas. This would have 
made the gas easier to strip ()Murakami fc Babul 119991) 
via ram pressure. Therefore, internal (SN), external (in- 
teractions, gas stripping, ram pressure) or global (reion- 
ization) mechanisms, can be at play, and the combina- 
tion of these might have shaped the morphology of some 
dSphs even b efore accretion into the satellite system of a 
large galaxy ([Sawala et al.l 12010). An analysis based on 
homogeneous SFHs derived for both the Milky Way and 
M31 satellites galaxies, together with reliable knowledge 
of their orbits, would be needed for a meaningful compar- 
ison of their properties with those of Tucana and Cetus, 
in order to gain some insight into the role of the environ- 
ment in the evolution of small systems. 

6.3. Tucana as a high-z galaxy 
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The range of estimated burst duration for Tucana 
(cr < 1.0 Gyr), allows us to infer the corresponding actual 
peak intensity of the initial burst: assuming cr > 0.1 Gyr, 
this ranges from 0.3 to 1.3 x 10~ 8 M Q pc~ 2 yr~ x . Note 
that this is the average value, as calculated over the whole 
area covered, from the center to ~ 3 core radii. Prelimi- 
nary analysis of the ip(t) as a function of radius (Hidalgo 
et al. in prep.) shows that the star formation rate in the 
innermost regions is « 4 times higher than the average 
value. This means that the SFR of Tucana, during the 
first stages of its evolution, was comparable to that ob- 
served for present da y Blue Compact Dwarf galaxies (see 
lAnnibali et al.ll2003l and references therein). 

Despite the fact that we can not put a stringent con- 
straint on the end of the star formation in Tucana, a solid 
result of our analysis is that Tucana was efficiently form- 
ing stars at ep ochs at least as old as z ~ 7. In various 
recen t works (iBouwens et al.l l2010al bl: iGonzalez et all 
I2010t lLabbe et all 120101; lOesch et all 120101 ), galaxies at 
similar redshift have been identified and characterized 
for the first time. However, all the galaxies so far de- 
tected have masses and luminos ity at least two order s of 
magnitude larger than Tucana ([Gonzalez et al.lf2010f ). 

7. SUMMARY AND CONCLUSIONS 

Deep HST/ACS data have allowed us to derive the life- 
time SFH of the Tucana dwarf galaxy, one of the most 
isolated dSphs in the LG. We have shown that Tucana 
experienced a strong event of star formation at the old- 
est possible age, > 12.5 Gyr ago. After the first ini- 
tial peak, the measured intensity of the star formation 
steadily decreased until stopping » 9 Gyr ago. The tests 
we performed with mock stellar populations disclose the 
broadening effect of the observational errors. We find 
that the actual underlying SFH is compatible with an 
episode of short duration, in the range cr < 1.0 Gyr if we 
assume a Gaussian profile ^(t) peaked 13 Gyr ago. Our 
attempt to put firm constraints on the age limits of the 
main event of star formation are hampered by the lim- 
ited time resolution, and thus we are not able to clearly 
answer the question of whether the reionization was de- 
cisive in ending Tucana's star formation activity. 

We explored alternative mechanisms (both external 
and internal) that may have shaped Tucana's evolution. 
On one hand, current measurements of its radial veloc- 
ity do not rule out the possibility that Tucana may have 
traversed the inner regions of the Local Group once. 
Therefore, a morphological transformation linked to a 
close interaction with a larger Local Group member at 
a time consistent with the end of its star formation can- 
not be ruled out. On the other hand, we explored the 
possibility that the feedback from supernova explosions 
might have been responsible for the gas loss. We used 
two different arguments to conclude that gas loss con- 
nected to supernova events must indeed have been very 
important in Tucana's early evolution: first, comparison 
of the total amount of energy released by the SNe ex- 
pected in the early evolution of Tucana with the model 
by Mc Low & Ferrara (1999) shows that Tucana could be 
in the blowout region. Second, investigating the chemical 



difference between two main sub-populations of distinct 
metallicity present in Tucana, we found evidence that the 
vast majority of the metals produced by the SNe must 
have been lost by Tucana to the inter-galactic medium. 

We devoted a particular effort to compare the proper- 
ties of Tucana and Cetus, the two isolated dSphs ana- 
lyzed consistently in the LCID project, that also share 
the important characteristic of being the most isolated 
dSphs in the LG. Both are (at first approximation) as old 
as the oldest Milky Way satellites, such as Draco, UMi or 
Sculptor, with no traces of star formation younger than ~ 
9 Gyr. The fact that they do not follow the morphology- 
density relation that has been observed in Milky Way 
dSph satellites poses interesting questions concerning the 
effectiveness of the environment in shaping the SFHs of 
dwarf galaxies. This, gives some support to models s uch 
as the one recently published by iSawala et all (|2010l ) in 
which internal mechanisms such as SNe, enhanced by 
the effects of cosmic reionization, are able to reproduce 
the main characteristics of dSph galaxies without having 
to invoke strong environmental effects. Still, new clues 
will come from a better understanding of the structure 
and kinematics of the stellar component of Tucana, and 
by comparing such properties in detail with those of the 
classical dSph satellites that were clearly affected by the 
environment. 

Despite the obvious similarities in the CMDs and SFHs 
of Cetus and Tucana, we also demonstrated important 
differences in their early evolution. We have shown that 
the formation time of the bulk of the stellar popula- 
tions in Cetus is clearly delayed compared to Tucana. 
This clearly appears in the derived SFHs, and other in- 
dependent indicators support the same conclusion: the 
morphology of the HB, the properties o f the RR Lyrae 
variable stars ([Bernard et all 1200 8. 2009) , and the char- 
acteristics of the RGB bump ([Monelli et alll2010al ). The 
most important conclusion we can draw from this com- 
parison is tha t it stro ngly reinforces the conclusions of 
I Monelli et all (|2010bl ). in particular that the vast ma- 
jority of the stars in Cetus were formed well after the 
end of the reionization epoch, therefore suggesting that 
the end of the star formation in Cetus was not predomi- 
nantly caused by it. This has important implications for 
state-of-the-art models on the effects of reionization in 
the early SFH of dwarf galaxies. 
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